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Motivated by the recent scanning tunneling microscopy (STM) and neutron scattering experi- 
ments, we investigate various charge density wave orders coexisting with superconductivity in the 
cuprate superconductors. The explicit expressions of the local density of states and its Fourier com- 
ponent at the ordering wavevector for the weak charge modulations are derived. It is shown that 
the STM experiments in BiiSr2CaCu20s+s cannot be explained by a site- or bond-centered charge 
modulation alone, but agree well with the presence of the dimerization hopping and transverse pair- 
ing modulations. We also calculate the spectral function for the charged stripes, which is measured 
by the ARPES experiments. 
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The coexistence of charged stripes and superconduc- 
tivity in high-T c superconductors has attracted a lot of 
both experimental and theoretical attention recently. In 
a novel STM experiment, Hoffman et al observed a four 
cell checkbroad local density of states (LDOS) modu- 
lation around the cores of superconducting vortices in 
Bi 2 Sr 2 CaCu20s+s by applying a magnetic field [1]. The 
charge modulation occurs in the Cu — O bond direction 
and in the energy range < E < \2meV. Very re- 
cently Howald et al discovered similar charge modulation 
at energy E = 2bmeV but in absence of magnetic field 
[2], The neutron scattering experiments on underdoped 
Y BCL2CU3O6.35 have also shown the charge density wave 
order with a period of eight lattice constants coexisting 
with superconductivity [3]. 

To explain the STM spectra in Bt2Sr2CaCu20s+8, a 
number of theoretical studies based on different models 
have been carried out by various authors [4-10]. How- 
ever, it is important to determine the correct charge den- 
sity wave orders in order to understand the origin of the 
charge stripes. In Ref. [8], Podolsky et al analysed in 
detail the influence of various patterns of translational 
symmetry breaking on the Fourier component of LDOS 
at the ordering wavevector Q by emyloping an approxi- 
mate technique. They concluded that the STM experi- 
ments [2] are consistent with the periodic modulation in 
the electron hopping. In contrast, by the numerical simu- 
lation of a d-wave superconductor with two-dimensional 
site charge density wave, bond charge density wave or 
pairing modulation, the pairing amplitude modulation 
comes closest to the experimental curves of Ref. [2] [9] . 
Therefore, there is no consensus about the constitution 
of the charged stripes. In this paper, in order to explain 
well the STM experiments, we solve strictly the d-wave 
superconductor with one-dimensional weak charge mod- 
ulations. The results show that the dimerization hop- 
ping and transverse pairing modulations are consistent 
with the STM experiments. We also discuss the ARPES 
experiments on the charged stripes, which verify the ex- 
istence of the hopping and pairing modulations. 

We start from the mean-field Hamiltonian of a d-wave 



superconductor 



#BCS = X] £ K c L C K<r + ^2 A K(<4 T C T _ Ki + C^KlCKT), 



Ker 



K 



(1) 



where £k = to + 1\ (cosK x + cosK y )/2 + t2CosK x cosK y + 
tz(cos2K x + cos2K y )/2 + t4,(cos2K x cosK y + 
cosK x cos2K y )/2 + t^cos2K x cos2K y and Ak = 
Ao(cos-Ka; — cosK y )/2. For nearly op- 

timally doped Bi2Sr2CaCu20s+s, io-5 = 
0.1305, -0.5951, 0.1636, -0.0519, -0.1117, 0.0510(eV) 
[11]. For underdoped YBa2Cu 3 Oe.35, io-2 = 
0.2445,-1.2000,0.3600 (eV) and t 3 _ 5 = 0.0(eV) (6 
percent doping). For both superconductors, we choose 
A = 0.0400ey. 

The charge modulations can be introduced phe- 
nomenologically into the Hamiltonian (1) by adding the 
charge ordered parameters, which have a general form 

H c = Ek ct (/kCk+q ct c Kct + /k c L c k+Q<t) 



Ek(.9k4c +qt 



5k c kt c -k-qi + h.c), 



(2) 

where Q = Qe^ (Q = ir/2 for S^S^CaC^Os+a 
and 7t/4 for Y Ba 2 Cu 3 6 .35), /k(= /-k-q) and gn{= 
,9-k-q) describe the hopping modulations and pair- 
ing modulations, respectively. We note that /k = 
Ai, A 2 e~ iQ / 2 , \ 3 cos(K x + Q/2)e" iQ / 2 and X 4 cosK y are 
the site-centered, bond-centered, longitudinal dimeriza- 
tion and transverse dimerization charge modulations 
while <7k = Agcos^a; + Q/2)e~ iC 3/ 2 and XecosK y are 
the longitudinal and transverse pairing modulations [8]. 
In this paper, we restrict our discussion to the case of 
weak charge modulations, i. e. small A's and do not con- 
sider the influence of the incommensurate spin orders to 
LDOS. In fact, the spin orders were not observed exper- 
imentally [2]. 

Taking the Bogoliubov transformation 



Ck+mQT 
L -k-mQJ. 



SkmO^kmO — CkmlV'kml 
£kml'0kmO + SkmO^kmlj 



(3) 



where k is restricted to the reduced Brillouin zone, 
m = 0,l,2,---,iV = (2tt/Q) - 1, = ±[1 + 

(-l)" e k+m.Q/-Ek+tn.Q], £kmo£krol = Ak+mQ^Sk+mQ 
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and E k+mQ = {el +mQ + A£ +mQ )2, the total Hamil- 
tonian H = Hbcs + Hq can be rewritten as 

H = Ekm I y( -1 )' y ^k+mQ'0k mI yV'km 1 / 
+ Ekrw{K+lm( k ) + /^r+lmW^L+l.VW + h.c}, 

where a^ +lm (k) = fk+ m Q[(-l) v+v Ckm+i^kmi/' - 

Ckm+ly+lCkmy' + l] and /3„ +lm (k) = g-k+mQ,[{-l) V ikrn+lv 
X £kroi/'+l + ( — ^-Y £km+li/+l£kmi/']- 

We define two-point Green's functions 

GZU^'-M = ~T< Tr[^ mv {r)A' m > v m >> (5) 

where T${t) denote the Fourier transform of 0(r) in 
Matsubara frequencies. Then the equations of motion 
for the Green's functions are 

[ioJn - ; (-l)^ k+mQ ]G^,(k, k';iw n ) 

- Y.v"(- l Y +U Ckmi^+ I /"|/k+(m-l)Q^m-l m '( k ; k '; iuJ n) 
+/k+mQ%+lm' ( k i k 'S 

_ E^C -1 )"" ^™+./"+i[sk+( m -i)Q%_V'( k . k '; 

+9k+mQCta'( k . k ';' W «)l = <W<W<W, 

(6) 

where 



^!,(k,k';; w „) = E,"(-i) (1 - I/> "am,+," 

xG:;(k,k' ;M „). 



(7) 



Obviously, combining Eqs.(6) and (7), the anomalous 
Green's functions in Eq. (5) can be solved by inverting 
a 2(N + 1) x 2(N + 1) matrix. To this end, we define 

£w( k ; iuj n) 



iu„-(-l)" E k+m q ' 

b m (k; iuj n ) = ^(-l)^km^km|,+iG° lly (k; iw n ), ^ 

C TO (k;iw n ) = Ei/ £kmiy+lG m „(k; JUn), 

^,(k,k';iw„) = (-l^-^'Ww&m^ 
x G^ w /(k;iw„), 

and introduce the 2(N + 1) vectors and g 



/ C 0l/ , 

/ y 0m' 



1 ,0° = 



g = 



,1 > 



\ ^ Nm' 
flu' 



UNm 



( ^Om' \ 



9om' \ 



9Nm' I 



(9) 



From Eqs. (6) and (7), wc obtain 

g = (i- M) _1 g, 

where / is 2(N + 1) x 2(7V + 1) unit matrix and 



M 



M 11 M 12 

M 21 M 22 



(10) 



(11) 



with M n ,M 12 ,M 21 and M 22 are (iV+1) x (iV+1) ma- 
trices whose matrix elements are as follows 

^mm' = ( a "i/k+mQ + ^flk+mQl^+l™' 

+ ( a m.fk+(m-l)Q + &mffk+( m — 1 m-lm' 5 

M mm i = { a m9k+mQ ~ ^/k+mQlWlm' 

+ ( a m.9k+(m-l)Q — ^m/k+(m-l)Q)^m-lro' j 
M mm i = (&m/k+ m Q + c m5k+ m Q)^m+lm' 

+ (&ro/k+(ro-l)Q + c m9k+(m-l)Q)^m-lm', 
M mm i = (^mffk+mQ — Cm/k+TnQ^m+lm' 

+ (^m.9k+(m-l)Q — c m/k+(m-l)Q)^m-lm' • 

(12) 

We note that the charged modulations observed in STM 
experiments are weak, i.e. /k+mQ and gu+mQ (or X's) 
are small. Expanding Eq. (10), we have Q = (I + M + 
M 2 + ■ ■ -)g. From Eq. (7), we finally obtain the green's 
functions 



c^(k,k'^„) = E,"(-i) (1 ^' > a m ,+^ 

x^,'(k,k';i W „). 



(13) 



To compare with the STM experiments, up to the first 
order in A's, we derive the local density of states 

p{r,u) = -^Im^h-T 7 < ^(^4^(0) >]L n _> a , +i0 + 

= -^f Im Ekm{ a m( k ;^«) + (/k + mQe ,Q,r + C.c) 

x [a m (k;iu} n )a m+1 (k;iuj n ) - 6 TO (k; iuj n )b m+1 (k; iu n )] 
+ (.9k+mQe iQ ' r + c.c)[a m (k;iuj n )b m+1 (k;iuj n ) 
+ &m(k; i<jj n )a m+ i (k; 

(14) 

where AT is the number of sites in the lattice and c r(J = 
^ 1/2 Ekm c k+mQa e l ( k+mC !) r . Obviously, the first term 
in Eq. (14) is nothing but LDOS for superconducting 
state. The other terms are the LDOS modulations with 
period 27r/Q due to the hopping and pairing modula- 
tions, respectively. The Fourier component of LDOS at 
the ordering wavevector Q is 

PQM - ^E r e- 4Q >(r,^) 

= "i;Ekm{/k+™q Im k( k ; iw„)a m+ i(k; iu> n ) 
-6 TO (k; iu n )b m+1 (k; iu) n )] 
+.9k+mQlm[a m (k; iu n )b m+1 (k; iu n ) 
+6 TO (k; iuj n )a m+ i(k; iw„)]}| iw?i _ >w+ ;o+ • 

(15) 

It is clear that when several of the charge density wave 
orders exist simultaneously, the total Pq(w) is obtained 
by a superposition of those of them. 

In Fig. 1, wc show the real parts of Pq(co) for dif- 
ferent charge density wave orders in BiiSriCaCuiO^+s, 
whose imaginary parts are zero or are proportional to its 
real parts. Obviously, the experimental curves (Fig. 3) of 
Rcf. [2] cannot be explained by the site- or bond-centered 
charge modulation alone (Fig. la). But our results for 
the dimerization hopping (Fig. lb, c) and transverse 
pairing modulations (Fig. le) are consistent with the 
STM experiments. Because the imaginary part of Pq(u)) 
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observed in the STM experiments is small, we conclude 
that the charge stripes are mainly formed by the trans- 
verse dimerization hopping and transverse pairing modu- 
lations. We note that our results are somewhat different 
from those of Ref. [8], this may be due to the approxi- 
mations using in Ref. [8]. We also note that the charge 
modulations discussed in Ref. [9] are two-dimensional 
but ours are one-dimensional. 

Fig. 2a-e show Repq(w) for YBa 2 Cu30e.35 with the 
same charge density wave orders as in Fig. 1. The curves 
in Fig. 1 and Fig. 2 are very different due to the different 
ordering wavevector and doping in two superconductors. 
For the dimerization hopping modulations (Fig. 2 b, c), 
Rc/5q(u;) has single peak rather than two peaks. For the 
pairing modulations (Fig. 2 d, e), there exist a strong 
peak and a weak peak at symmetric positions. These 
results are expected to be verified by the STM experi- 
ments. 

Finally, we discuss the ARPES experiments on the 
charge stripes, which measure the spectral function. 
From Eqs. (3) and (13), we obtain the spectral function 
up to the second order correction 

A k (lu) = A m (k,uj) 

= -±Im[-T< c k +mQT(-r)c k+roQT (0) >]L„^ +40 + 
= -ilma m (k;io;„){l + X; y [am+ ! ,-i(k;ia;„)/* +(m+i/ _ 1)Q 
+&m+^-i(k; «w„)5 k+ ^ m+iy _ 1 ^g][a m+I/ (k; iu n )fk+( m +is-i)Q 
+b rn+l ,(k; iuJ n )9k+(m+v-i)Q] + [&m-i(k; iUn)f H+( m -i)Q 
+c m _i(k; «w„)g£ +(m _ 1)Q ][a m (k; iWn)#k+(m-i)Q 

-6 TO (k; iw n )/k+(m-i)Q] + [&m+i(k; iu n )fa+mQ 
+c TO+ i(k; iw„).g k+mQ ] [a m (k; «w„)5 k+mQ 
— b m (k; *w„)/ k+m q]}| iw?i ^ w+i o+ 
= A k(uj) + 5A k (lj) 

(16) 

The first term in the above equation is the spectral func- 
tion of the superconducting state. The others are those 
due to the weak charge modulations, which are the sec- 
ond order small quantities. Obviously, the spectral func- 
tions Ak(w) for different hopping (pairing) modulations 
have similar energy dependence at the same momentum 
K. So it is difficult to clarify which kinds of hopping 
(pairing) modulations the charged stripes belong by the 
ARPES experiments. However, the hopping and pairing 
modulations have unique energy dependence for Ak(w) 
which could be distinguished experimentally. Fig. 3 and 
Fig. 4 show the spectral functions Ak(w) at various 
momenta K for Bi2Sr2CaCu20s+s and F-BCI2CU3O6.35 
with hopping and pairing modulations, respectively. 

In summary, we study the effects of various charge den- 
sity wave orders on the d-wave superconductor. We con- 
clude that the dimerization hopping and transverse pair- 
ing modulations and experimental results are in excellent 
agreement. The origin of such charge modulations is due 
to the hopping and gap disorders in the superconductor, 
which will be investigated in detail in another paper [12]. 



We also discuss the ARPES experiments on the charged 
stripes, which can distinguish the hopping and pairing 
modulations. 
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FIG. 1. Energy dependence RepQ(a;) of the Fourier com- 
ponent of the LDOS at Q = (tt/2,0) for Bi 2 Sr 2 CaCu 2 8+s 
with different charge density wave orders, (a) A site- or 
bond-centered charge density wave, (b) A longitudinal dimer- 
ization charge density wave, (c) A transverse dimerization 
charge density wave, (d) A longitudinal pairing modulation, 
(c) A transverse pairing modulation. In both (b) and (d), 
Ropq(lo) was multiplied by -1 to compare conveniently with 
the experiments. 

FIG. 2. Repq(w) at Q = (tt/4,0) for Y Ba 2 Cu i 6 . 35 with 
the same charge density wave orders as in Fig. 1. 

FIG. 3. Spectral function Ak(ui) vs uj at various K for 
Bi 2 Sr 2 CaCu 2 Og + s with hopping and pairing modulations. 

FIG. 4. Spectral function Ak(w) vs to at various K for 
Y Ba 2 CuzO§.zh with hopping and pairing modulations. 
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